Deflagration with quantum and dipolar effects in a model of a molecular magnet 
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Combination of the thermal effet in magnetic deflagration with resonance spin tunneling controlled 
by the dipole-dipole interaction in molecular magnets leads to the increase of the deflagration speed 
in the dipolar window near tunneling resonances. 
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Molecular magnets such as Mni2 Ac, possessing an 
effective large spin S^IO, are famous as mesoscopic 
systems demonstrating magnetic bistability due to the 
strong uniaxial anisotropy. Spectacular spin tunneling 
under the anisotropy barrier in molecular magnets was 
first seen in the steps in dynamical hysteresis curves,- 
The steps correspond to the values of the longitudinal 
magnetic field B z at which the energy levels in both po- 
tential wells match. Here resonance spin tunneling leads 
to a faster relaxation responsible for a step of the magne- 
tization. To the contrast, off resonance the main channel 
of relaxation is thermal activation over the top of the 
barrier. The difference between the two types of relax- 
ation is shown in Fig. [TJ In fact, spin tunneling requires 
a transverse field or any other term in the Hamiltonian 
that breaks the axial symmetry. Pure spin tunneling in 
the right panel in Fig. [1] requires that these terms be 
sufficiently strong, such as the transverse field of about 3 
T in Mni2. In the case of weaker tunneling interactions, 
the intermediate situation of a thermally assisted tunnel- 
ing is realized. In this case spins tunnel after thermally 
mounting up to below the top of the barrier The role 
of tunneling in the case of weaker tunneling interactions 
can be interpreted as some lowering of the barrier near 
resonances. For stronger tunneling interactions, the bar- 
rier is removed completely at resonances. 

Tunneling and relaxation in molecular magnets can be 
described by the density matrix equation^ the most com- 
prehensive account of which is given in Ref. @. The 
latter numerically implements the universal spin-phonon 
interaction suggested in Refs. 0H- This interaction is 
due to distortionless rotation of the crystal field acting 
on the spins by transverse phonons and it is completely 
expressed in terms of the crystal-field Hamiltonian Ha 
without any unknown spin-lattice coupling constants. 

Experiments of 2005 by Myriam Sarachik group 
showed the existence of propagating deflagration (burn- 
ing) fronts in the molecular magnet Mni2 Ac that are 
similar to chemical burning. 9,10 Javier Tejada group ob- 
served peaks in the deflagration speed on the bias mag- 
netic field B z that were interpreted as contribution of 
resonance spin tunneling ^ A detailed, mainly classical, 
theory of the magnetic deflagration including the ignition 
threshold and the accurate prefactor in the Arrhenius- 
type expression for the speed of the burning front was 
proposed in Ref. [ll|. 

The physics of deflagration is based on triggering re- 



laxation from a metastable state over potential barrier 
by the temperature increase as the result of relaxation 
accompanied by energy release. The burning front forms 
because the temperature in the regions still unburned 
(e.g., before the front) rises as the result of heat conduc- 
tion from the hot areas where burning just occured. The 
two main ingredients of deflagration thus are the Arrhe- 
nius dependence of the relaxation rate on temperature 
(making burning in the cold areas before the front negli- 
gibly slow) and heat conduction. Deflagration is mathe- 
matically described by the system of coupled i) rate equa- 
tion for the number of particles (magnetic molecules) in 
the metastable state and ii) the heat conduction equa- 
tion. 




FIG. 1: Thermal activation in the nonresonant case vs tun- 
neling in the resonant case. 

Subsequent theoretical quest for an essentially quan- 
tum form of deflagration lead to the discovery of self- 
organized fronts of tunneling, a non-thermal process trig- 
gered by the dipolar field (rather than by temperature) 
that can bring the system on or off resonance.— ~— A 
hallmark of these fronts is the self-consistent adjustment 
of the metastable population (or magnetization) to the 
optimal spatial profile that creates the dipolar field that 
is constant in some region of space and brings the system 
on resonance. The width of the resonance region form- 
ing the front core is about the transverse dimension R 
of the sample that allows an efficient tunneling and thus 
front propagation. On the other hand, before and after 
the front core the system is off resonance and tunneling 
is blocked. Fronts of tunneling can be realized in the 
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dipolar window of the external field B z 
< B z - B k < B[ D) . 



(1) 



Here B k is the field corresponding to the kth resonance 
in the absence of the dipolar field and Bz is the dipo- 
lar field created by the uniformly magnetized elongated 
crystal. It was showr>12. that the adjustment mechanism 
is robust with respect to resonance spread (e.g., due to 
defects) smaller than B Z D \ 

The aim of the present paper is to unify the theories of 
the standard (hot) deflagration 2 ^ and fronts of tunneling 
(cold or quantum deflagration) 1 13 i 14 

In the sequel we will use the generic model of a molec- 
ular magnet with Ha = ~DS Z , where the tunneling res- 
onance fields are given by 



B k = kD/( gf i B ), fc = 0,l,.. 



(2) 



Resonance tunneling occurs at B to t.z = B z + B^ ps B k 
between the metastable ground state — S) and an excited 
state at the other side of the barier \m!) with m! — S — k. 
At temperatures much smaller than the barrier height 
(e.g. at the temperature of the deflagration front) one 
can describe magnetic molecules as two-level systems oc- 
cupying the states \±S) . Let us denote the probability 
for a molecule to be in the metastable state |— S) as n. 
Then the average value of the effective pseudospin o~ z is 



a z = 1 - 2n, 



(3) 



so that n = 1 corresponds to a z = —I. The general 
expression for the longitudinal component of the dipolar 
field on magnetic molecule i is the sum over positions of 
all other moleculs j 
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Here vq is the unit-cell volume, D z 
field, and 



is the reduced dipolar 
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Inside a uniformly magnetized ellipsoid, a z = const, the 
dipolar field is uniform and one has D zz = D zz a z , where 



D zz = Dt h) + 4^ (l/3 



(0) 



v is the number of magnetic molecules per unit cell (v = 2 
for Mni2 Ac) and = 0, 1/3, and 1 for a cylinder, 
sphere, and disc, respectively. The reduced dipolar field 
in a sphere D z s z h ^ depends on the lattice structure. For 
Mni2 Ac lattice summation yields D^z = 2.155 that 
results in D z c z yl) = 10.53 for a cylinder. Then Eq. g]) 
yields the dipolar field B [ Z D) ~ 52.6 mT in an elongated 
sample that was also obtained experimentally^ 

For simplicty we consider a long crystal of cylindrical 
shape of length L and radius R with the symmetry axis 
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FIG. 2: Relaxation rate T(B Z ,T) in a generic model of a 
molecular magnet. 



z along the easy axis, magnetized with o~ z = a z {z). The 
latter assumption makes the problem tractable numeri- 
cally. In this case the reduced magnetic field along the 
symmetry axis has the form 1 ^— 



D zz {z) = 



dz'- 



2-KvR 2 a z (z') 

I i3/2 

(z 1 - zf + R 2 



ka z (z), (7) 



where 



3tti73 - D<f^ = 4irv - > 0, 



(8) 



k = 14.6 for Mn 12 Ac. For other shapes such as 
elongated rectangular, one obtains qualitatively similar 
expressions Now the total field is given by 
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z {z) = B z +Bi D \z) = B 



■D zz (z). (9) 



One of the dynamical equations of the model is the 
relaxation equation for the metastable population n(t, z) 



dn(t, z) 
dt 



-T(B totiZ (z),T(z)) \n(t,z)-n^(T) 



(10) 

In Eq. (TTU1) r (. . .) is the relaxation rate taking into ac- 
count both thermal actication over the barrier and reso- 
nance spin tunneling that is calculated from the density 
matrix equation. As i?tot,z depends on n(z) everywhere 
in the sample via Eqs. (|7|) and ©, this is an integro- 
differential equation. In the sequel we will set => 
that is a good approximation for strong enough bias. 2 

The second equation is the heat conduction equation 
that is convenient to write in terms of the energy £ of 
the system per unit cell as in Ref. In the full-burning 
case n( cq ) = this equation has the form 



d£(t, , 
dt 



d_ d£(t,z) 
dz dz 



n AE 



dn(t 7 z) 
dt 



(11) 



In Eq. (ITT|) k is the thermal diffusivity and AE is 
the energy released in the transition of a spin from the 
metastable state to the ground state, 



AE = 4hDS 2 , 



2DS 



(12) 
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The relation between the energy £ and temperature is 
given by £{T) = / Q T C(T')dT', where C(T) is the exper- 
imentally measured heat capacity per unit cellJ£ 

To solve the system of Eqs. (ITU1) and (jlll) numerically, 
it is convenient to introduce reduced variables^ 



£ = 



£ 



n AE 



tT f , 



r 

r = F' 

Id 



(13) 



where no < 1 is the initial population of the metastable 
state and Tf is the relaxation rate at the flame temper- 
ature Tf defined by the energy balance n AE = £(Tf) 
and some fixed value of B to t.z that we set to the reso- 
nance field Bk- The characteristic distance Id = WKf/Tf 
defines the width of the deflagration front in the case of 
normal (thermal) deflagration and Kt is the thermal dif- 
fusivity at Tf . In terms of these variables, Eqs. (ITU1) and 
(fTTjl become 



d£ 

Ih 
On 
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d _ d£ dn 
dz dz dr 

V ( B tot , z ,T(£) 



(14) 
(15) 



where f = T/Tf is the reduced relaxation rate and 
K = k/ Kf. It remains to add the expression for B to t, z 
in reduced variables, Eq. (|9]) with D zz (z) given by Eq. 
([7]) with z z and R => R = R/ld- The important 
parameter R is the ratio of the width of the front of tun- 
neling that is of order R (see Refs. lil^flil ) to the width 
of the standard deflagration front Id^- 

Eqs. (fl~4")l and ([15)) are solved numerically by choos- 
ing a finite-length sample and discretizing the problem 
in z. This yields a system of ordinary differential equa- 
tions in time. We set k — 1 for simplicity. Before solving 
the equations, T was calculated from the density ma- 
trix equation^ for the transverse field B± = 3 T and 
tabulated as a function of B to t,z and £ . As T increases 
by many orders of magnitude near tunneling resonances, 
one has to use many different values of B to t.z for inter- 
polation here. In Fig. [5] one can see that for such a 
strong transverse field the barrier is reduced to zero at 
resonance where T practically does not depend of tem- 
perature. Thus near the resonance the cold deflagration 
should dominate, while off resonance the regular defla- 
gration should take place. 

For the discussion it is convenient to consider the en- 
ergy bias W — £-s — £ m i between the two resonant levels, 



W = (S + m 1 ) gn B (B a + B^ - B fc ) = W c , 
It is convenient to use the reduced external bias 
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W {D) . 
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where Ed — (Sg/J-s) /^o is the dipolar energy, Eu/Ub = 
0.0671 K for Mn 12 Ac. At the right end of the dipolar 
window of Eq. (JTJ> one has B z = B k +B^ . Thus with the 
help of Eq. gj one obtains W cyA = (1/2) (1 + m'/S) D zz , 




FIG. 3: Profiles of the metastable population n and the total 
bias W across the front for two values of the external bias 
Wext: (a) M4xt = 0, laminar regime; (b) Woxt = 5, non- 
laminar regime 



i.e., Wext ~ Dzz for small bias, m! « S. We will see that 
in the case of strong tunneling the speed of the quantum 
deflagration front has a maximum at the right end of the 
dipolar window, W ext ~ D zz yV> = 10.53 for Mni 2 Ac. 

Cold deflagration can be ignited by the field sweep 
across the resonance. In this case ignition occurs around 
the ,, magic" value W ext — 5 that corresponds to B z — 
Bk = 22 mT i 13 ' 14 Outside the dipolar window fronts of 
tunneling do not exist. On the other hand, standard de- 
flagration can be initiated, at any bias, by a quick tem- 
perature rise on one side of the sampled Applying this 
metnod of ignition here, we will see that within the dipo- 
lar window the process is modified by spin tunneling and 
the speed of the burning front can significantly increase 
for R > 1, especially at the right end of the window. 

There are two regimes of propagation of non-thermal 
fronts of tunneling: Laminar and non-laminar. Lami- 
nar regime with a smooth front takes place in the left 
part of the dipolar window, < B z — B k < 10 mT (or 
< Woxt < 1-3), while the non-laminar regime with 
frozen-in quasiperiodic spatial patterns of the magneti- 
zation behind the front is realized in the right part of 
the dipolar window. In both regimes burning is not com- 
plete and becomes less complete with increasing the bias. 
In the laminar regime the residual magnetization and 



4 



k=l,B ± = 3T 



Normal 
deflagration 




Dipolar window 



K = 1(1 



Normal 
deflagration 



N 52 mT W„ 



FIG. 4: Reduced front speed v vs bias field for different values 
of the reduced transverse size R = R/ld, Id being the width 
of the thermal deflagration front. For such a strong applied 
transverse field, the effect of tunneling is dramatic. 



the front speed were calculated analytically. 14 The front 
speed increases with the bias. In the non-laminar regime, 
quasi-periodic frozen-in patterns of magnetization dete- 
riorate the resonance condition, and the front speed de- 
creases with the bias after the breakdown of the laminar 
regime (see Fig. 5 of Ref. \l4 ). Thermal mechanism of 
deflagration leads to complete burning of this residual 
metastable population that smoothens the dipolar field 
profile in the sample and improves the resonance condi- 
tion inside the front core. This leads to the increase of 
the front speed because of spin tunneling in the whole 
dipolar window. 

Results of numerical calculations for the spatial pro- 
files of the metastable population n and the total bias 



W in the front for R = 1 arc shown in Fig. [3J Both 
in the laminar and non-laminar regimes, there is a re- 
gion where W = and resonant tunneling takes place, 
causing a greater slope of n(z). Behind the frond (on the 
left) metastable population n burns to zero via thermal 
mechanism. 

Numerical solutions for the reduced front speed v = 
v/(ldTf) (Ref. [ll|) for the generic model with B± = 3 
T and R = 1 and 10 are shown in Fig. [4] Within the 
dipolar window the front speed can largely exceed the 
speed of regular deflagration and depends on the trans- 
verse crystal size R parametrized by R. At B± = 3 T 
the maximal values of v are attained at the right end 
of the dipolar window, followed by a steep drop towards 
the standard-deflagration result outside the dipolar win- 
dow. For smaller transverse fields such as 2 T, the effect 
of spin tunneling is weaker and v reaches a maximum 
somewhere in the middle of the dipolar window, depend- 
ing on R. 

Measurements of the speed of deflagration fronts*^ 
were done in zero or small transverse field, so that the 
influence of resonance spin tunneling on the front speed 
is not so dramatic as in Fig. 2) It would be highly in- 
teresting to perform deflagration experiments in strong 
enough transverse field to see the big effect of tunneling 
on the front propagation. Changing thermal contact of 
the crystal with the environment, one can boost or sup- 
press the thermal mechanism of magnetic burning thus 
isolating thermal and quantum effects from each other. 

Numerous useful discussions with E. M. Chudnovsky 
are greatfully acknowledged. 

This work has been supported by the NSF Grant No. 
DMR-0703639. 



1 J. R. Friedman, M. P. Sarachik, J. Tejada, and R. Ziolo, 
Phys. Rev. Lett. 76, 3830 (1996). 

2 D. A. Garanin and E. M. Chudnovsky, Phys. Rev. B 56, 
11102 (1997). 

3 E. M. Chudnovsky and D. A. Garanin, Phys. Rev. Lett. 
79, 4469 (1997). 

4 D. A. Garanin, X. M. Hidalgo, and E. M. Chudnovsky, 
Phys. Rev. B 57, 13639 (1998). 

5 D. A. Garanin and E. M. Chudnovsky, Phys. Rev. B 59, 
3671 (1999). 

6 D. A. Garanin, arXiv:0805.0391 (2008). 

7 E. M. Chudnovsky, Phys. Rev. Lett. 92, 120405 (2004). 

8 E. M. Chudnovsky, D. A. Garanin, and R. Schilling, Phys. 
Rev. B 72, 094426 (2005). 

9 Y. Suzuki, M. P. Sarachik, E. M. Chudnovsky, S. McHugh, 
R. Gonzalez-Rubio, N. Avraham, Y. Myasoedov, E. Zel- 
dov, H. Shtrikman, N. E. Chakov and G. Christou, Phys. 



Rev. Lett. 95, 147201 (2005). 

10 A. Hernandez-Minguez, J. M. Hernandez, F. Macia, A. 
Garcia-Santiago, J. Tejada, P. V. Santos, Phys. Rev. Lett. 
95, 217205 (2005). 

11 D. A. Garanin and E. M. Chudnovsky, Phys. Rev. B 76, 
054410 (2007). 

12 D. A. Garanin and E. M. Chudnovsky, Phys. Rev. B 78, 
174425 (2008). 

13 D. A. Garanin and E. M. Chudnovsky, Phys. Rev. Lett. 
102, 097206 (2009). 

14 D. A. Garanin, Phys. Rev. B 80, 014406 (2009). 

15 S. McHugh, R. Jaafar, M. P. Sarachik, Y. Myasoedov, H. 
Shtrikman, E. Zeldov, R. Bagai, and G. Christou, Phys. 
Rev. B 79, 052404 (2009). 

16 A. M. Gomes, M. A. Novak, R. Sessoli, A. Caneschi, and 
D. Gatteschi, Phys. Rev. B 57, 5021 (1998). 



